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We study the critical behavior of the random g-state Potts model in the large-q limit on the 
diamond hierarchical lattice with an effective dimensionality d c ff > 2. By varying the temperature 
and the strength of the frustration the system has a phase transition line between the paramagnetic 
and the ferromagnetic phases which is controlled by four different fixed points. According to our 
renormalization group study the phase-boundary in the vicinity of the multicritical point is self- 
similar, it is well represented by a logarithmic spiral. We expect infinite number of reentrances in 
the thermodynamic limit, consequently one can not define standard thermodynamic phases in this 
region. 



I. INTRODUCTION 

Disorder is an inevitable feature of real materials and it 
could result in dramatic changes in the physical proper- 
ties of the systems, in particular in the vicinity of phase- 
transition points. The effect of disorder can be even more 
pronounced when both random ferromagnetic and anti- 
ferromagnetic couplings are present. Then the combined 
effect of frustration as well as thermal and disorder fluc- 
tuations could lead to a rich phase diagram, which con- 
tains a spin-glass (SG) phase and different critical and 
multicritical points^. 

Among the models studied in this context the random 
bond q-state Potts model 2 - - — plays an important role. For 
q = 2 we have the well known Edwards- Anderson model^ 
of Ising SG and for q = 3 and q = 4 it is used to describe 
orientational and quadropular glasses^, respectively. For 
large values of q the model is suggested to be a plausible 
model of supercooled liquids. This relation seems to be 
correct in the infinite-range model^, but for a system 
with finite-range interactions further studies are needed 
to clarify this point. On the cubic lattice for q = 3 n i 12 
and later for q = 5 and 6^ numerical evidence was ob- 
tained in favor of a thermodynamic transition into a sta- 
ble SG phase. For q = 10 on the cubic lattice in the 
first studies no SG phase has been observe d 1 ^ 12 . Re- 
peating the calculation on a one-dimensional chain with 
long-range power-law interactions, which is equivalent to 
a short-ranged system below the upper critical dimen- 
sion, evidence is obtained for a SG phase at low enough 
temperature^. Continuing this study further for even 
larger q-s seems to be very difficult, in particular due to 
the very narrow range of temperature where the presence 
of the SG phase is expected. 

In this paper, in order to gain some insight about 
the possible behavior of the large-g model we use the 
Migdal-Kadanoff renormalization group (RG) method^, 
which provides exact results on hierarchical lattices^. 
This method has already been used for the Potts 
model (with finite value of q) for random ferromagnetic 



couplings^—, while for random fcrro- and antiferromag- 
netic couplings this method is combined with the Gaus- 
sian approximation^!. Studies of the model in the large-q 
limit has been performed by one of the present authors 
in Rcfi^ 2 - to which we will refer as paper I. In paper I the 
lattice with an effective dimensionality, d c g = 2 is stud- 
ied in more detail, in which case the critical properties of 
the model are calculated by solving numerically a set of 
integral equations. 



In this paper we continue these studies when the dia- 
mond hierarchical lattice has a larger branching number, 
b > 2, which corresponds to d e g > 2. In this case the 
RG equations are more complicated, which are then iter- 
ated numerically. For different values of b we determine 
the phase diagram in the temperature-frustration plane, 
which is found to contain a ferromagnetic and a para- 
magnetic phase, but there is no SG phase. The phase 
transition is shown to be controlled by a random ferro- 
magnet (RF) fixed point when the driving force of the 
transition is the temperature, whereas it is controlled 
by a zero temperature (Z) fixed point when the driving 
force is frustration. We calculate the critical exponents 
in the two attractive fixed points for different effective 
dimensionalities. The attractive regimes of the RF and 
Z fixed points in the phase transition line is separated by 
an unstable multicritical (MC) fixed point, the vicinity 
of which is studied in more detail. Our numerical inves- 
tigations predict that the phase-diagram around MC has 
a singular structure. 



The structure of the paper is the following. The model 
and the method of investigation is described in Secllll 
Results about the phase diagram and the critical prop- 
erties of the system are presented in Sec lIIIi while the 
vicinity of the MC point is studied in Sec lIVl Our paper 
is closed by a Conclusion in the final section. 
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II. MODEL AND MIGDAL-KADANOFF 
RENORMALIZATION 

We consider the g-state Potts model denned by the 
Hamiltonian: 



The initial condition is given by 



Inq 



(1) 



in terms of the Potts spin- variables Sj = 1,2, ...,q at 
site i of the lattice. Here <5 Si . s is the Kronecker-symbol, 
the summation runs over nearest-neighbor pairs and the 
couplings Jij are independent and identically distributed 
random numbers, which can be either positive or nega- 
tive. The lng in the denominator insures that the tran- 
sition temperature stays finite when q goes to infinity^. 
In the random cluster representation of the model the 
partition function Z is dominated by one diagram: 



subleading terms 



(2) 



and the free energy is given by F = —<f>/0, where /3 = 1/T 
with the convention ks = 1. 

In this paper the Potts model is placed on the diamond 
hierarchical lattice, which is constructed recursively from 
a single link. At each step a link is replaced by a unit, 
which consists of b parallel branches, each branch con- 
taining two bonds in series. At generation n, the length 
L n measured by the number of bonds between the two 
extreme sites A and B is L n = 2", and the total number 
of bonds is 



B n = (2b) n = Z#*W with d eS {b) = 



ln(26) 
ln2 



(3) 



where d c ff(b) is an effective dimensionality. 

In the following we consider fixed-spin boundary con- 
ditions: when the two extreme sites A and B are in the 
same (different) state the partition function is denoted 
by Z*' 1 (Z^ 2 ). Their ratio is given by: 



-2-1,2 



(4) 



where I n = /3F™ tcr and F}^ 1 ' = F 7 \> 2 - F 1 / 1 is the in- 
terface free energy. The scaled interface free energy I n 
satisfies the recursion cquatior 



,17.19 



? w = n 



q^+^+^-2) 



(5) 



which for large-g becomes g-independent (sec paper I) : 



7 . , = \^ $ ll^ I^- 

J-n+l — / j ~x y- n , ^ n 

i=l 

Here the auxiliary function is: 



(0) 



7« , = max r 7 (i) + 7(2) ; i) _ max Ui) 5 7(2) } ^ 



r(i) _ 



(8) 



where Ji is the value of the ith coupling. 



III. PHASE DIAGRAM AND CRITICAL 
PROPERTIES 

A. Numerical pool method 

The phase diagram of the random system is studied 
numerically for different numbers of the branching num- 
ber b > 2, which corresponds to an effective dimension 
(i c ft > 2. The original distribution of the couplings is 
given in a box-like form: 



V{J) 



if j£— < J < -r^—, 

otherwise. 



(9) 



(7) 



with p < 1 as already used in Paper I for 6 = 2. For p > 
all couplings are random ferromagnetic and in the limit 
p — > 1 we have the pure system. For p < there are also 
negative bonds, their fraction is increasing with decreas- 
ing p. In the numerical calculations we have used the so- 
called pool method. Starting with N random variables 
taken from the original distribution in Eq. (|9]), we gen- 
erate a new set of N variables through renormalization 
at a fixed temperature, T, using Eqs. ^ and J7]). These 
are the elements of the pool at the first generation, which 
are then used as input for the next renormalization step. 
We check the properties of the pool at each renormaliza- 
tion by calculating the distribution of the scaled interface 
free energy, its average and its variance. In practice we 
have used a pool of N = 5 x 10 6 elements and we went 
up to n ^ 70 — 80 iterations. We mention that in order 
to study the sensitivity of the results on the form of the 
initial distribution we have also used Gaussian- and sym- 
metric bimodal distributions instead of Eq.©, but the 
structure of the phase diagram and the properties of the 
fixed points remain the same. 



B. Phases 

At a given point of the phase diagram, (p,T), the 
rcnormalizcd parameters display two different behaviors, 
which are governed by two trivial fixed points. 

i) In the paramagnetic phase the scaled interface free 
energy renormalizes to zero, thus its average is I = and 
its variance is A7 = 0. 

ii) In the ferromagnetic phase the scaled interface free 
energy goes to infinity. Its average behaves at a length, L 
as I ~ L ds , with d s = d c ff — 1 being the effective dimen- 
sion of the interface. Also the variance of I is divergent: 
AI ~ L , with a droplet exponent 8 > 0. The length L 
has been defined by Eq[3J 
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We note that in the large-g limit there is no spin-glass 
phase, even at zero-temperature, contrary to the known 
results for q = 2 and q = 3 . 24 ' 25 . 




FIG. 1: Phase-diagrams for different branching numbers: 
6 = 2, 3, 4, 8 and in the b — ► oo limit. The system is in 
the paramagnetic phase for high temperature (large T and r) 
and/or for large frustration (small p), otherwise it is in the 
ferromagnetic phase. The phase transition lines are indicated 
by red and green, along which the transition is controlled by 
the RF and the Z fixed points, respectively. The red and 
green lines meet at the MC points where there are several 
re-entrances for b > 2. 



0.75 




0.25 



FIG. 2: (Color online) Probability distribution of the scaled 
interface free energy at the RF fixed point for b = 2, b = 3, 
6 = 4 and 6 = 8. At / = there is an extra delta peak with 
respective strength 0.12808, 0.25052, 0.31889, and 0.43609. 



the fixed point is given in Figj2]for 6 = 2, 3, 4, and 8. 
At this fixed point the distribution contains no negative 
/ values, but there is a 6-peak at I = 0, the strength 
of which is increasing with 6, see the caption of FigOJ 
In the vicinity of the transition line at a distance AT = 
T c — T in the ferromagnetic phase the average value and 
the variance of / behave as: 



The two phases are separated by a phase transition 
line T c (p) which can be calculated accurately for a given 
pool and its true value can be obtained by averaging over 
different pools and taking the N — >• oo limit. The phase 
diagrams are shown in Fig. Q] for different values of the 
branching number in the plane r = (1 +p)T/J versus 
p. (In the vertical axis we used a nonsingular scaling 
combination, since the average coupling: J = 2±i_ is zero 
at p = — 1.) For the case b oo, the exact expressions of 
the transition lines can be obtained and these are given in 
Appendix. Interestingly the phase diagram shows several 
reentrances in the regime p < 0. (see also Fig. [5]). This 
behavior is different from that observed at 6 = 2 (in which 
case only one reentrance is present) and we are going to 
study its properties in more detail in Sec. IIVI 



I(AT,L) 
A/(AT, L) 



L 



Cav(AT) 

L 



Cvar(AT) 



+ 



(10) 



Here the correlation lengths are divergent as: £ av (AT) ~ 
(AT)-"" and £ V ar(AT) - (AT)""™, respectively. 

iii) The other part of the transition line between MC 
and Z (indicated by green line in Fig. [T]) is controlled 
by the Z fixed point, which is located at T = t = 0. 
Here the interface free-energy grows with the size as: 
I{L) = L 6z uj, where the ui are 0(1) random numbers 
and the droplet exponent is 6z > 0. Here we use the 
scaled variable, i n = I n /I n , which obeys the renormal- 
ization group equation from Eq.([5]): 



C. Fixed points 

For all values of b the phase-transition line is found to 
be controlled by four different fixed points. 

i) The fixed point of the pure system (P) is located at 
p = 1 and t = 2(26 — l)/6 and it describes a first-order 
transition. This fixed point is unstable for any small 
amount of disorder, see paper I. 

ii) Along the transition line for p < 1 between P and 
the multicritical point MC (indicated by red line in Fig. 
Q]) the transition is controlled by the RF fixed point. Here 
7 = 0(1) and AI = 0(1) and the distribution P(I) at 



I, 



+ i 



In 



with a n+ i = I n+ i/I n and 



(11) 



(i^+iW o)-inax(i« *< a >,0) 



(12) 

The fixed-point value of the ratio: a z is related to the 
droplet exponent as: 0z = log a z / log 2. 

The distribution of the scaled variable: i is shown in 
Fig|3]for various values of 6. In the vicinity of the transi- 
tion line at a distance AT = T c — T in the ferromagnetic 
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FIG. 3: (Color online) Probability distribution of the scaled 
interface free energy at the Z fixed point. At i — there is 
a delta peak with strength 0.1728 10' 
for 6 = 3, and 0.940010" 5 for 6 = 4 
pool of size 5 x 10 6 we have found no zero value. 



3 for b = 2, 0.492810 -4 
For 6 = 8 and with a 



phase the average value and the variance of the interface 
free energy scales as: 



I (AT, L) = 
AI(AT, L) = 



L" 



[£ W (AT)]«<.-** 



Kvar(AT)] 



(13) 



respectively. Here the correlation lengths, £ av (AT) and 
£ V ar(AT) are divergent at AT = 0, with the critical ex- 
ponents: i/f v and z^ r , respectively. 

iv) The attractive regions of the RF and Z fixed points 
are separated by the unstable MC fixed point. Here the 
interface free energy behaves as: I = 0(1) and AI = 
0(1), similarly to the RF fixed point. The distribution, 
P(I) at the MC fixed point is shown in FigU]for various 
values of b. This distribution has a finite jump at / = 
for b > 2, which has interesting consequences about the 
structure of the phase diagram at MC, what we are going 
to analyze in SeclIVI 



D. Critical exponents 

The critical exponents for b > 2 are calculated numeri- 
cally by the pool method. In the calculation of the RF ex- 
ponents we have first fixed the parameter of the disorder 
distribution, p, then calculated the location of the critical 
point: T c (p). For a given pool T c (p) has been calculated 
accurately, with a precision of 1CP 12 . Then estimates for 
the critical exponents are calculated from the relations in 
Eqs. fTU]) . At the Z fixed point the temperature is T = 
therefore we have to calculate its position: pz- Then the 
critical exponents are calculated from (IT31 . Repeating 
the procedures for different parameters and for different 
pools we have obtained estimates of the set of critical 



FIG. 4: (Color online) Probability distribution of the scaled 
interface free energy at the MC fixed point for 6 = 2, 3, 4 
and 8. For 6 > 2 there is a discontinuity at / = (a step 
representation of the probability distribution has been used 
to stress this feature). In addition there is a delta peak at 
1 = with respective weight 0.0057, 0.012854, 0.014713 and 
0.016626 for 6 = 3, 6 = 4 and 6 = 8. 



exponents with a relative precision of about 10~ 3 . We 
have observed that the correlation length critical expo- 
nents calculated from the average and from the variance 
are identical for a given fixed point, i.e. v av = v vaT = v 



and similarly 



This relation is found 



to be exact before for b = 2 (sec paper I) and now it is 
demonstrated numerically for other values of the branch- 
ing number. We note that for finite value of q the two 
correlation length exponents are generally different, at 
least on the same diamond hierarchical lattice^. The set 
of numerically calculated critical exponents are collected 
in Table U for different values of the branching number. 
For completeness we also present here the numerically 
exact results obtained for b = 2 in paper I. 



TABLE I: Critical exponents of the disordered Potts model in 
the large-g limit on the diamond hierarchical lattice. Results 
for branching number 6 = 2 are from paper I. 
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v z 


1.729 


1.258 


1.051 


0.738 


0.555 


0.440 



IV. MULTICRITICAL POINT 

The phase-diagram in the vicinity of the MC point is 
different for 6 = 2 and for b > 2, see in Fig IT] In the 
former case there is just one reentrance of the phase- 
boundary and at the MC fixed point there is traditional 
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multicritical behavior, see in paper I. On the contrary for 
b > 2 the phase diagram around MC has more than one 
reentrance. We have explored the phase diagram around 
MC by using larger and larger magnifications and the re- 
sults for b = 8 are presented in FigJS] Interestingly the 
phase-diagram shows a self-similar structure: the reen- 
trances are repeated in different scales. Let us denote 
the location of the MC point by: r* and p* and then 
use relative positions of the phase-boundary, which are 
given by Ap = p c — p* and At = t c — r* . According 
to the results in Fig(5] by enlarging the scale by roughly 
a factor of ten the phase-diagram has approximately the 
same reentrant structure. With the pool method, having 
a large, but finite number of realizations we could go up 
to a magnification of po « 2 x 10 -7 in the phase diagram, 
up to which scale one could observe at least seven similar 
reentrant structures. 









FIG. 5: Phase-diagram in the vicinity of the MP point: the 
framed regions are enlarged in the next figure. Note the (self- 
similarity of the phase boundary in different scales. The 
green (bright) region represents the paramagnetic phase, the 
red (dark) region is the ferromagnetic phase. 



It is instructive to use polar coordinates, p and </>, 
and plot the phase-diagram in a semi-logarithmic scale 
as a function of ln(p/ p$) and (f>, which is presented in 
FigH] In this figure the paramagnetic and the ferromag- 
netic phases are represented by two almost perfect spi- 



rals, which means that the phase-boundary around the 
MC point is approximately given by a logarithmic spiral. 
We note that for other values of the branching number 
(b = 3 and 6 = 4) we have observed a similar structure 
of the phase diagram around the MC point. 




FIG. 6: Phase-diagram in the vicinity of the MP point in 
logarithmic polar coordinates, see text. 



We expect that the size of the elementary magnifica- 
tion step, po, is decreasing when the size of the pool is 
increasing and in the thermodynamic limit the magnifi- 
cation step goes to zero. In this case the phase diagram 
at the MC point becomes singular. Passing through the 
MC point the phase diagram has infinite number of reen- 
trances and sufficiently close to MC one can not define 
a stable thermodynamic phase, since its stability range 
measured either in temperature or in concentration of 
couplings goes to zero. 



V. CONCLUSION 

In this paper we considered the g-state Potts model 
for large q with randomly mixed ferro- and antiferromag- 
netic couplings and studied its critical properties on the 
diamond hierarchical lattice with an effective dimension- 
ality, d e g > 2. By using the Migdal-Kadanoff renormal- 
ization scheme, which is exact for this lattice, the phase 
diagram is explored numerically using the pool method. 
As for d c f[ = 2 the phase diagram contains a ferromag- 
netic and a paramagnetic phase, but there is no spin glass 
phase, even at zero temperature. The phase transition in 
the random system is controlled cither by a random ferro- 
magnct fixed point, in which the driving force is the tem- 
perature, or by a zero temperature fixed point, in which 
the transition is due to frustration. The two regimes are 
separated by a multicritical point in the vicinity of which 
the system shows unusual properties. 

Both in the RF and the Z fixed points the transition is 
continuous and we have calculated numerically the crit- 
ical exponents, which are collected in Table [J] The cor- 
relation length exponents calculated from the scaling of 
the average or the variance are found to be equal. At the 
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MC point the distribution of the reduced interface free 
energy is a discontinuous function, having a finite jump 
at I = 0. At the same time the phase-diagram close to 
MC has a self-similar, logarithmic spiral-like structure. 
Consequently the phase boundary itself has a singularity 
and in the traditional sense one can not define a stable 
thermodynamic phase here. 

The studies presented in this paper can be extended 
in different directions. First, it would be interesting to 
study by analytical methods the behavior of the system 
in the vicinity of the MC point for b > 2. i) To prove or 
disprove, that the elementary magnification step at the 
MC point, po> approaches zero in the thermodynamic 
limit, ii) To show rigorously, that the distribution func- 
tion, P(I), in the MC point is strictly discontinuous, as 
shown in FigfJ] iii) Also would be interesting to clarify 
the possible relation between these two expected proper- 
ties. 

Another extention of the present study is to investi- 
gate, if the singular behavior at the MC point could ex- 
ist in another systems, too. As far as the random Potts 
model on hypcrcubic lattices is concerned it has the same 
fixed point structure for finite value of q, thus generally 
there is an MC point in these model o 24 ' 26 . Probably in 
sufficiently large dimension (d = 3 could be enough), and 
for large-g, the MC singularity could exist. Numerically, 
however, it would be very difficult to explore the details 
of the process, as already seen in the problem of detec- 
tion of the SG phas o 13 ' 14 . A discontinuous distribution 
of some observables at this point, such as the reduced in- 
terface free energy, could be the sign of such singularity. 
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Appendix: The phase-diagram for b — >• do 

The phase-diagram for b — > oo can be analytically 
treated, here we consider the initial distribution in 
Eq.©. First we note that in this limit in the recursion 



equation Eq. §6§ there is a sum of a large number of ran- 
dom variables and the renormalized value of the interface 
free energy depends on the sign of the average value of 

the auxiliary function: /^(/q 1 ^, Iq 2 ^)) = ($)• Here we 

can differentiate between three cases. 

i) If (<f>) > then all the renormalized parameters are 
I± ^> 1, consequently we are at the F fixed point and 
due to the central limit theorem the droplet exponent is 
9 = d s /2. 

ii) If (<I>) = 0, then all the 1\ values are equal to zero, 
thus the P fixed point has been reached. 

iii) Finally, if ($) < 0, then all the renormalized pa- 
rameters are I\ <C 0, thus in the next renormalization 
step we have I? = and again the P fixed point has been 
reached. 

In the next step we compute the sign of ($) in the 
different parts of the phase diagram in FigJT] Let us 
denote the end-points of the box-distribution in Eq.© 
by x = -£ and y = ^. In the region y < 1 one has 
($) = and only for x < and y > 1 we can have 
($) < 0. A direct calculation gives in this region ($) 
as the ratio of a third degree polynomial in x and y by 
(y — x) 2 and the numerator has two different expressions 
according to the sign of (x+y+1). Then expressing x and 
y with p and r, one gets that the transition lines shown 
in Fig. [T] The upper boundary is a horizontal line with 
r = 4. The lower boundary is given by two expressions: 

-3r 3 + 12t 2 - 24r + 32 + 12pr 2 - 48pr + 48p = 
-Up 3 - 48p 2 + 24p 2 r - t 3 + 16 = , 

which are valid for p < p and p > p, respectively. The 
two lines join at: 

p = i (1 - 4cos6>) 

f = ^(l-cos0) 

arctan (3\/7) — tt 
9 = 3 ' 

and the extremities of the lower boundary line are at 
p = — oo, r = 2 and p — 2 cos (|tt) — 1, t = 0. 

We note that in the limit b — > oo the ferromagnetic 
phase between 2 < r < 4 extends to all values of p. For 
large, but finite value of 6, however, there is a paramag- 
netic phase for p < — \/b. 



* Electronic address: igloi.ferenc@wigner.mta.hu 
1 For reviews on spin glasses, see K. Binder and A. P. Young, 
Rev. Mod. Phys. 58, 801 (1986); M. Mezard, G. Parisi, and 
M. A. Virasoro, Spin Glass Theory and Beyond (World Sci- 
entific, Singapore, 1987); K. H. Fisher and J. A. Hertz, 
Spin Glasses (Cambridge University Press, Cambridge, 



1991); Spin Glasses and Random Fields, edited by A. P. 
Young (World Scientific, Singapore, 1998). 

2 D. Elderfield and D. Sherrington, J. Phys. C 16, L497 
(1983). 

3 D. J. Gross, I. Kanter, and H. Sompolinsky, Phys. Rev. 
Lett. 55, 304 (1985). 



7 



4 H.-O. Carmesin and K. Binder, J. Phys. A 21, 4053 (1988). 

5 M. Scheucher, J. D. Reger, K. Binder, and A. P. Young, 
Phys. Rev. B 42, 6881 (1990). 

6 G. Schreider and J. D. Reger, J. Phys. A 28, 317 (1995). 

7 O. Dillmann, W. Janke, and K. Binder, J. Stat. Phys. 92, 
57 (1998). 

8 S. F. Edwards and P.W. Anderson, J. Phys. F: Met. Phys. 
5, 965 (1975). 

9 K. Binder and J. D. Reger, Adv. Phys. 41, 547 (1992). 

10 G. Parisi, in Complex Behavior of Glassy Systems ed M 
Rubi and C Perez- Vicente (Berlin, Springer, 1997) p 111; 
S. Franz and G. Parisi Physica A 261, 317 (1998); M. 
Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999); 
M. Mezard Physica A 265 352 (1999); G. Parisi, Physica 
A 280 115 (2000); X. Xia and P. G. Wolynes Phys. Rev. 
Lett. 86, 5526 (2001). 

11 C. Brangian, W. Kob, and K. Binder, J. Phys. A 36, 10847 
(2003). 

12 L. W. Lee, H. G. Katzgraber, and A. P. Young, Phys. Rev. 
B 74, 104416 (2006). 

13 R. A. Banos et al., J. Stat. Mech. P05002 (2010). 

14 R. S. Andrist, D. Larson, and H. G. Katzgraber, Phys. 
Rev. E 83, 030106(R) (2011). 

15 A. A. Migdal, Sov. Phys. JETP 42, 743 1976; L. P. 
Kadanoff, Ann. Phys. N.Y. 100, 359 (1976). 

16 A. N. Berker and S. Ostlund, J. Phys. C 12, 4961 (1979); 



M. Kaufman and R. B. Griffiths, Phys. Rev. B 24, 496 

(1981) ; R. B. Griffiths and M. Kaufman, ibid. 26, 5022 

(1982) . 

W. Kinzel and E. Domany, Phys. Rev. B 23, 3421 (1981). 

B. Derrida and E. Gardner, J. Phys. A 17, 3223 (1984); B. 
Derrida, in Critical Phenomena, Random Systems, Gauge 
Theories, edited by K. Osterwalder and R. Stora, Pro- 
ceedings of the Les Houches Summer School of Theoretical 
Physics, 1984 (North- Holland, Amsterdam, 1986), p. 989. 

D. Andelman and A. N. Berker, Phys. Rev. B 29, 2630 
(1984). 

C. Monthus and T. Garel, Phys. Rev. B 77, 134416 (2008). 
A. Benyoussef and M. Loulidi, Phys. Rev. B 53, 8215 
(1996). 

F. Igloi and L. Turban, Phys. Rev. B 80, 134201 (2009). 
J.-Ch. Angles d'Auriac and F. Igloi, Phys. Rev. Lett. 90, 
190601 (2003); M.-T. Mercaldo, J.-Ch. Angles d'Auriac, 
and F. Igloi, Phys. Rev. E 73, 026126 (2006). 

E. S. Sorensen, M. J. P. Gingras, and D. A. Huse, Euro- 
phys. Lett. 44, 504 (1998). 

I. Morgenstern and K. Binder, Phys. Rev. B22, 288 (1980), 
I. Morgenstern and H. Horner, Phys. Rev. B 25, 504 
(1982), W. L. McMillan, Phys. Rev. B28, 5216 (1983). 
J.L. Jacobsen, and M. Picco, Phys. Rev. E 65, 026113 
(2002). 



